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Abstract We have examined the structure of hot accretion flow with a large-scale mag- 
netic field. The importance of outflow/wind and thermal conduction on the self-similar 
structure of a hot accretion flows has been investigated. In comparison to the accretion 
disk without winds/outflow, our results show that the radial and rotational velocities of 
the disk become faster however it become cooler because of the angular momentum and 
energy flux which are taking away by the winds/outflows, but thermal conduction opposes 
the effect of winds/outflows not only decrease the rotational velocity but also increase the 
radial velocity as well as the sound speed of the disk. In addition we have studied the ef- 
fect of global magnetic field on the structure of the disk. We have found out that all three 
components of magnetic field have a noticeable effect on the structure of the disk such as 
velocities and vertical thickness of the disk. 
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1 INTRODUCTION 

Black hole accretion disks provide the most powerful energy -production mechanism in the universe. It 
is well accepted that many astrophysical systems are powered by black hole accretion. Most of the ob- 
servational features of black hole accretion systems can be explained through The standard thin disks ( 
geometrically thin and optically thick accretion disks) with a significant success (Shakura & Syunyaev, 
1973). The ultraviolet and optical emission which is observed in quasars frequently has a role in the 
thermal radiation in the standard disks (SD) that surrounds by the massive black holes in quasars (e.g., 
Sun & Malkan, 1989). On the other hands, the SD seems unable to recreate the spectral energy distribu- 
tions (SED) of in a lot of other luminance sources (e.g., Sgr A.) accreting in low rates, and the advection 
dominated accretion flows (ADAF) were recommended to be present in these sources (Narayan & Yi, 
1994, 1995). In the ADAF model, the most energy which is released from the infalling gases in the 
accretion flow is changed to the internal energy of the gas. only a small fraction of the energy in ADAFs 
gets radiated away, thus, their radiation efficiency is a lot less than that for SD (see Narayan et al., 1998, 
for a review and references therein, Kato, Fukue & Mineshige 2008). 

Outflows/wind have been found in numerical simulations of hot accretion flow by both magnetohy- 
drodynamic (MHD) and hydrodynamic (HD) simulations. In some pioneer paper outflow was reporte- 
din HD simulations are by Igumenshchev & Abramowicz (1999, 2000) and Stone, Pringle & Begelman 
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(1999). Igumenshchev & Abramowicz (2000) have shown that convective accretion flows and flows 
with large-scale circulations have significant outward-directed energy fluxes, which have important im- 
plications for the spectra and luminosities of accreting black holes. Most recent work focusing on the 
convective outflows is by Yuan & Bu (2010). As a result, they have shown that the mass accretion rate 
decreases inward, i.e., only a small fraction of accretion gas can fall onto the black hole, while the rest 
circulates in the convective eddies or lost in convective outflows. Stone & Pringle (2001) in their MHD 
simulations have reported the appearance of wind and outflow. In their simulations the net mass accre- 
tion rate is small compared to the mass inflow and outflow rates at large radii associated with turbulent 
eddies. 

There are a lots of models indicate that modeling the hot accretion flows is a big challenging and 
it is controversial problem which should be investigated. Thermal conduction is one of the important 
neglected physical phenomena in the modeling of ADAFs. Recent observations of the hot accretion flow 
around AGNs indicated that they should be based on collision-less regime (Tanaka & Menou 2006). 
Chandra observations provide tight constraints on the physical parameters of gas. Tanaka & Menou 
(2006) have used these constraints ( Loewenstein et al. 2001 ; Baganoff et al. 2003; Di Matteo et al. 2003; 
Ho et al. 2003) to calculate the mean free path for the gas materials. They have reported and suggested 
that thermal conduction can be as a possible choice as a mechanism for transporting sufficiently extra 
heating is provided in hot advection dominated accretion flows. So it should be important to consider 
the role of thermal conduction in a ADAF solution. Shadmehri (2008), Abbassi et al. 2008,2010, Tanaka 
& Menou (2006), Faghei 2012, have studied the effect of hot accretion flow with thermal conduction 
with a semi-analytical method; dynamics of such systems have been studied in simulation models (e. g. 
Sharma et al. 2008; Wu et al. 2010). Shadmehri (2008) has shown the thermal conduction oppose the 
rotational velocity, but increase the temperature. Abbassi et al. (2008) have shown that for this problems 
there are two types of solution, high and low accretion rate. They had plotted the radial velocity for both 
of the solutions which have shown that it will modified by thermal conduction. 

previous studies about black hole accretion disks recommend that a large-scale magnetic field root 
in the ISM or even central engine, will be taken inward and compressed closed in the black hole by the 
accreting materials (Bisnovatyi-Kogan & Ruzmaikin 1974, 1976). Therefor, large-scale magnetic field 
has a great role on the structure of hot flow because of highly ionized flow. The effect of magnetic field 
on the structure of ADAFs were also reported by ( Balbus & hawley 1998, Kaburiki 2000, Shadmehri 
2004, Meier 2005, Shadmehri & khajenabi 2005, 2006, Ghanbari et al. 2007, Abbassi et al. 2008, 2009, 
Bu,Yuan & Xie (2009)). Akizuki & fukue 2006 and Abbassi et al 2008, presented the self-similar solu- 
tions of the flow based on the vertically integrated equations. At the same time, they have stressed on 
the intermediate case where the magnetic force is compared with other forces by this assumption that 
the physical variables of an accretion disk have just radius function and also they explained more about 
disk toroidal magnetic field . 

In this paper, we first extend the work of Akizuki & Fukue 2006; Zhang & Dai 2008 and Abbassi 
et al. 2008, 2010 by considering a general large-scale magnetic field in all the three components in 
cylindrical coordinates (r, ip, z) and then discuss effects of the global magnetic field on the flows with 
thermal conduction and wind. We adopt the treatment that the flow variables are functions of the disk 
radius, neglect the different structure in the vertical direction except for the z-component momentum 
equation. We also compare our results with pervious studies, in which a large-scale magnetic field, 
thermal conduction and wind are neglected. 

2 BASIC EQUATIONS 

We are interested in analyzing the structure of a magnetized ADAF with a large-scale magnetic field 
where thermal conduction and wind play an important role in energy and angular momentum trans- 
portation. Then, we consider a accreting and rotating disk around a compact black hole of mass M*. 
So, for a steady axi-symmetric accretion flow,(i.e., d/dt = d/dip = 0), we write our equations in 
the cylindrical coordinates (r, ip, z). We vertically integrated the equations and, then all our physical 
variables become only function of radial distance r. Moreover the disk is supposed to have Newtonian 
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gravity in radial direction and also we neglect the general relativistic effects, we suppose that the disc 
have magnetic field with three components B r , B v and B z . 
The continuity equation reads, 

where v r is radial velocity and S is the surface density at the cylindrical radius r, which we can define 
as £ = 2pH, H is the half-thickness and p being midplane density of the disk. The mass loss rate by 
wind/outflow is represented by M w , so 

M w (r) = I 4itr'rh w (r')dr' (2) 



Here m w (r) is the mass-loss rate per unit area from each disk face. According to Blandford & Begelman 
1999; Shadmehri 2008, We write the dependence of accretion rate as follows 

M = -2irrEv r = A/ ( — ) s (3) 
ro 

where Mq is the mass accretion rate of the disk material at the outer edge of the disk r$ and s with s are 
constant of order unity (Blandford & Begelman 1999). For a disk with outflow/wind, s > while in the 
absence of wind/outflow, we consider s = (Fukue 2004). 
By considering equations (d}-©, we can have 

fn w = sp,{Ly-> (4 ) 

For the equation of motion in the radial direction we can write, 

Here, c s the isothermal sound speed, which it can define as c? s = p gas /p, Pgas being the gas pressure and 
Dp is the rotational velocity. Also, c r , c v and c z are Alfven sound speeds in three direction of cylindrical 
coordinate and define as, z = j (47rp). 

We suppose that, just (^-component of viscose stress tensor is important which is define as t r(p = 
firdil/dr , where p(= vp) is the viscosity and v is the kinematic constant of viscosity. So we have 

v = ac s H (6) 

Here a is the standard viscous parameter which is less than unity and constant (Shakura & Sunyaev 
1973). 

The angular transfer equation with considering outflow/wind can be written as 

?L*-(rv ) = — — (i^Er 3 — ) + — (VSc ) + ^Ss. - ^ d ^ w (7) 
r dr v r 2 S dr dr y 7 ^ dr v r 27rS dr 

The last term on the right side of above equation represents angular momentum carried by the 
wind/outflow. As we know, I = corresponds to a non rotating outflow/wind and I = 1 to outflow- 
ing material that carries the specific angular momentum away(e.g., Knigge 1999). 
Similarity, by integrating over z of hydrostatic balance , we will have 

nlH* - -±= Cr ± { ^ Cz)H = cl + \(cl + 4) (8) 
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Now we write the energy equation considering heating and cooling processes in an hot accretion 
flow. We suppose that the generation energy due to viscous dissipation and are balanced by the advection 
cooling and energy loss of wind/outflow. Thus 

2Hv r cfJL = fullr^f - 2 4UrF s ) - W)<M (9) 



B v = 




v r B v ), 


B z = 


d , 
- T (v r B z ) 
dr 


v r B z 
r 



7 — 1 dr s dr dr r dr 2 

where the second term on right hand side of the energy equation represents transfer of energy due to the 
thermal conduction and according to Cowie & Makee 1977, F s (= b<f> s pc 3 s ) will be the saturated conduc- 
tion flux. The coefficient <f> a is less than unity. Here, the last term on right hand side of energy equation 
shows the energy loss due to outflow/wind. In our case, we consider the dimensionless parameter 77 it a 
free parameter and constant (Knigge 1999). 

To measure the magnetic field escaping/creating rate, we can write the three components of induc- 
tion equation, (B r , B v , B z ), as 

B r = 0, (10) 

(ID 
(12) 

Here B TtiPtZ is the field creating or scaping rate due to dynamo effect or magnetic instability. Now 
we have a set of magnetohydrodynamical equations (MHD), that describe the structure of magnetized 
advection dominated accretion flows (ADAFs). The solutions of these equations are related to viscosity, 
magnetic field strength (3 r ,ip,z, the advection parameter /, the thermal conduction parameter <fi s and 
wind/outflow parameter s of the disks. In the next section, we will seek a self-similar solution for our 
main equations. 

3 SELF-SIMILAR SOLUTIONS 

To better understanding of the physical processes taking place in our hot accretion disk, we seek self- 
similar solutions of the above equations. Therefore, following Narayan & Yi (1994), we write similarity 
solutions as 

E(r) = coE ( — ) a "* (13) 
ro 



GM*.r 1 

v r (r) = -c 1 J ( — ) ^ (14) 



= e^/—(-) > (15) 

c^(r) = c 3 (^)(^)- 1 (16) 
B 2 r ^ z GM* 

— _ Zp r ip Z C3( 

47T/9 r r 



= = 2/3,,^,,c 3 (^)(-)- 1 (17) 



H(r)=c 4 r (-) (18) 

Here ,cq, c\, C2, C3 and C4 are constants and will be determined later, ro and So are made the equations 
in non-dimensional forms. Moreover the constants /3 r>tp>z measure the radio of the magnetic pressure in 
three direction to the pressure of gas (i.e., /3 r . v . z = p m ag,r,<p,z/Pgas)- 
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The field creating/scaping rate Br.^.z is assume to be as 

Br,tp,z = B r Q^ 0z0 ( — ) 2y 2 

r 



(19) 



where -B,o,v>o,zO are constant. 

If we substitute the above self-similar solutions in to the main disk equation such as, the continuity, 
momentum, angular momentum, hydrostatic balance and energy equation of the disk, we can obtain the 
following system of equations to be solved for cq, c\,C2 c 3 and C4: 



1 - [( f 



CoCi 
3 



2 ) + (*-f)A + (* + ^)/V|c3 



1. 



2 ClC2 



Sl — C2 
CO 



C4= 2 



(S - ^) 2 PrP Z c\ + 4(1 + p r + /^)c 3 + (S - -)VMzC 3 



1 



where m is the non dimensional mass accretion rate which can write as 



9. 2; — _ , , § s m 

-afc 2 ^/c 3 c 4 - 5(j) s (s - 2)c 3 - - -77— 

4 4 C 



(20) 

(21) 
(22) 

(23) 
(24) 



M 



27rr So\/GM*/ro 
Also, for the magnetic field creating/scaping rate, we have, 

Bar = 



7 GM* AttcqCs'So 



f'4 



^02 



1. 3, GM» /47r/3 z coc 3 So 
2( s -2 )Cl 77TV ~ 4 



(25) 

(26) 
(27) 

(28) 



Without mass outflow/wind, thermal conduction and magnetic field,(i.e. s = l = r) = <f) s = f3 r = fi l p = 
/3 Z = 0), our equations and their solutions are reduced to the standard ADAFs solution (Narayan & Yi 
1994). Moreover in the absence of outflow/wind and thermal conduction they are reduced to Zhang & 
Dai 2008. 

Now we can analysis behavior of the solutions in the presence of large-scale magnetic field, out- 
flow/wind and thermal conduction. The parameters of this model are the advection parameter /, the 
radio of the specific heats 7, the standard viscose parameter a, the mass loss parameter s, the degree of 
magnetic pressure to the pressure of gas , f} r , /3 V and (3 Z , the conduction parameter (j> s and I, rj parame- 
ters corresponding to the non-rotation and rotation wind/outflow and the loss of energy by wind/outflow. 



Four panels in figure[T]show the variations of coefficients cq, c\, c 2 and c 3 in term of wind parameter 
s for different values of thermal conduction parameter <fr s , i. e., (j> s = 0.001 (dotted line), (f> s = 0.01 
(dashed line) and <f) s = 0.02 (solid line), corresponding to a = 0.2, 7 = 4/3, /3 r = f3 v = j3 z = 1.0, 
77 = I = 0.1 and / = 1.0 (fully advection). In upper-left panel of figure[T]we can see that for non-zero s, 
the surface density is lower than the standard ADAF solution and for stronger outflows, the reduction of 
the surface density is more evident . Also surface density decreases by increasing of thermal conduction. 
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Fig. 1 Numerical coefficient c;s as a function of wind parameter s for several values of 
thermal conduction parameter (f> s . The dotted, dashed and solid lines correspond to (f> s = 
0.001, 0.01 and 0.02 respectively. Parameters are set asa = 0.2, 7 = 4/3, j3 r = fi v = /3 Z = 
1.0, 77 = I = 0.1 and / = 1. 



Upper-right panel of figure Q] represents the behavior of radial infall velocity to the Keplerian one. 
Although radial velocity is slower than the Keplerian velocity, it becomes faster by increasing / and a. 
We know that the value of s measure the strength of outflow/wind, and large values of s denotes strong 
wind/outflow. So we can see that, the redial flows of the accretion materials become faster by increasing 
of wind/outflow parameter s. Moreover thermal conduction increases the accretion velocity. It is clear 
from bottom panels of FigureQ] the rotational velocity increase by adding wind parameter, s. Although 
thermal conduction decreases the rotational speed. It is shown that temperature decreases for strong 
outflows (down-right panel). On the other hand, outflow/wind plays as a cooling agent (Shadmehri 
2008; Abbassi etal2010). 

Figure|2]shows how the coefficients c;s depend on the magnetic parameter in radial direction j3 r for 
different values of outflow/wind parameter s. We can see that the surface density become larger if the 
wind parameter s(0 — 0.3) increases. It is obvious that s = represent no-wind solutions (Tanaka & 
Menou 2006). By adding f3 r , we see that the surface density and sound speed decrease. While, the radial 
and rotational velocity increase. The same as figure Q] outflow increase radial and rotational velocity, 
although the sound speed and surface density decreases by increasing of wind parameter. 

In figure|3]the behavior of the coefficients as versus toroidal magnetic field are plotted for different 
values of wind parameter s. From the upper-left panels of figure|3] we can see that for P v {= — 5), the 
accretion velocity is sub-Keplerian, and that becomes faster by increasing of wind/outflow parameter s. 
Moreover, when the toroidal magnetic field becomes stronger, the accretion velocity of flow increases. 
This is because the magnetic tension terms, which dominates the magnetic pressure term in the radial 
momentum equation that assists the motion of accreting materials in radial direction. The down-right 
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Fig. 2 Numerical coefficient c,s as a function of magnetic parameter (j r for several values 
of wind parameter s. The dotted, dashed, dash-dotted and solid lines correspond to s = 
0.0, 0.1, 0.2 and 0.3 respectively. Parameters are set as a = 0.2, 7 = 4/3, /3 V — /3 Z = 1.0, 
tj> s = 0.01, Tf = 1 = 0.1 and / = 1. 



panel in figure|3]display the rotational velocity of accretion disc. We can see that for the given advection 
parameter /(= 1.0) and also j3 r = f3 z = 1.0, the rotational speed increase as the toroidal magnetic 
field become larger. The considerable matter about this panel is that, when the toroidal magnetic field is 
very strong, the rotational velocity will be Super-Keplerian (i.e., /3 V ~ 5). Moreover the surface density 
of the disk decreases when the toroidal magnetic field become large. In figure [3] the isothermal sound 
speed of the disk has distinctive properties between f} v < 1 and f3 v > 1. When fi v is below unity, the 
isothermal sound speed decrease, although for j3 v > 0, it increases. It means that the is a minimum 
near (3 V ~ 1, because the parameters f3 r = /3 Z (= 1.0) are present in our calculation and both of them 
decrease the isothermal sound speed. So when the toroidal magnetic filed dominate, sound speed start 
be become large. 

In figure|4] we have shown that the coefficients as with vertical component of magnetic field for 
several values of s. We figured out that a strong (/3 Z ) leads to an decrease of the accretion velocity v r , 
rotational speed v v as well as isothermal sound speed c s of the disc. That is because of a large magnetic 
pressure in the z direction prevents the accretion flow from being accreted, and will decreases the effect 
of gas pressure as accretion proceeds. Although the surface density of the disk increases by increasing 
of ^-component of magnetic field. 

4 GENERAL PROPERTIES OF ACCRETION FLOWS 

In previous section it has introduced self-similar solutions and their properties for a advection dominated 
accretion flow bathed in general large-scale magnetic field with thermal conduction and outflow/wind. 
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Fig. 3 Numerical coefficient as as a function 
of s. The dotted, dashed, dash-dotted and solid 
respectively. Parameters are set as a = 0.2, 7 = 
and / = 1. 
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of magnetic parameter j3 v for several values 
lines correspond to s = 0.0, 0.1, 0.2 and 0.3 
4/3, P r = p z = 1.0, (j) s = 0.01, rj = I = 0.1 



Having this self-similar solutions allow us to investigate the general properties of the disk, so we can 
see how thermal conduction and wind/outflow will effect on the physical quantities of the disks. 

The mass-accretion rate in the disk can be written by a simple lapse function of radius due to the 
effect of outflow and thermal conduction. Using the self-similar solution we can have the mass accretion 
rate as, 

M = -2nrZv r = 2^c ciS (^^)5 (-) s = M out ( — y (29) 

ro r r out 

where r out and M out are the outer's radius and accretion rate respectively. As we have in last section 
thermal conduction, wind and large scale B-field will effect on cq, c\ through their parameters s, <j> and 
Pr,<p,z- In the case of an accretion disk with no wind, s = 0, the accretion rate is independent of the 
radial distance, while for those with wind/outflow, s > 0, the accretion rate decreases with radius as 
we expect. Because winds start from different radii, the mass-loss rate is not fixed and it depends on 
the radius. Then, some region of accretion flow are not concentrated at the center, but are expand over a 
wide space. 

If the central star is a black hole of mass M„, it will grows via accretion. It simply can estimate the 
BH growth rate using accretion rate which also will affected by wind/outflow, thermal conduction and 
large-scale B-field, implicitly. 

We can see the radial behavior of temperature of the present self-similar advection dominated ac- 
cretion discs with a general large-scale magnetic field, thermal conduction and outflow. Optical depths 
of a hot accretion flow are highly dependent on their accretion rate. These ADAFs occur in two regimes 
which depend on optical depth and their mass accretion rate. In a high mass accretion rate, the optical 
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Fig. 4 Numerical coefficient qs as a function of magnetic parameter (3 Z for several values 
of s. The dotted, dashed, dash-dotted and solid lines correspond to s = 0.0, 0.1, 0.2 and 0.3 
respectively. Parameters are set as a = 0.2, 7 = 4/3, (i r — f3 v = 1.0, <j> s = 0.01, rj = l = 0.1 
and / = 1. 



depth becomes very high and the radiation generated by accretion flow can be trapped into the disc. 
This type of accretion disc is known as slim disc or optically thick disc. In the low mass-accretion rates, 
the disc becomes optically thin and the cooling time of accretion flow is longer than the accretion time- 
scale. So, the generation energy by accretion flows remains mostly in the discs and the discs are not able 
to radiate their energy efficiently. Then, we should expect that the Large-scale magnetic field will effect 
on the optical depth and accretion rate implicitly through the coefficient c u , C\. 

We try to show that the radiative appearance of optically thin discs, (such as L v ~) , is not able to be 
calculated in ease. Also it should be demonstrated the significance of advective cooling in optically thin 
discs. To clarify the matter, we have to state how cooling and heating change with /, r and M. In all, 
the radiative cooling was so complicate in case of dependency parameters. In optical thin discs, the gas 
emission does not have a black body continuum. Spectra emission seems to be ruled by Bremsstrahlung 
cooling by non relativistic or relativistic, synchrotron and also Compton cooling, to deal will the opti- 
cally thick case, where we are trying to study it, the radiation pressure dominates and sound speed is 
related to radiation pressure. Then, we write the average flux F as: 

F = oTt = £-n gas = 3 {E ^&^y-i (30) 
8H 8 C4 Tq ro 

Here, H gas = Sc^ is the height-integrated gas pressure, T c is the disc central temperature and a is 
constant of the Stefan-Boltzmann . in order to calculate the disk optical thickness in vertical direction 
of the optically thick case we can write: 

t = tUs = is oCo (-) s -^ (31) 
2 2 r 
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where k is the electron-scattering opacity. Then, we calculate the effective flux of the disc surface as: 

<rT* ff = ^ = (^L E )ir-i (32) 
1 ' t Wwa C4 

Here, Le is the Eddington luminosity of the central object which is define as L E = 47rc^^. If we 
radially integrate these equations, we can have the disc luminosity as, 

L disc = [° Ut 2aT* f 2nrdr = h E ~H — ) (33) 

we can not emphasize more that the effective temperature and the luminosity of the disc are not in- 
fluenced by the magnetic field, thermal conduction and either mass-loss through outflow (there are 
no ^ r ,4>,z, </>s and s dependencies). But wind/outflow might affect the radiative appearance of the disc 
through the C3 and C4 in these formulae, implicitly. When we have mass-loss wind/outflow compared 
with the case of no mass-loss, then the average flux decline through the disc . For the mass-loss case, 
the surface density and therefore optical depth decrease and we can see the deep inside of the disc. 



5 CONCLUSIONS 

In this paper, we studied an accretion disk in the advection dominated regime by considering a large- 
scale magnetic field and in presence of wind/outflow and thermal conduction . Also, some approxima- 
tions were made to simplify our main equations. We supposed an static and axially symmetric disc with 
the a-prescription of viscosity, v = ac s H and a set of similarity solution was presented for such a 
configuration. We have extended Akizuki & Fukue 2006; Zhang & Dai 2008 and Abbassi et al 2010 self 
similar solutions to present the structure of the advection dominated accretion flows (ADAFs). Also, we 
ignored self-gravity of the discs and the general relativistic effects. 

Our results have shown that strong wind/outflow can have lower temperature and there are satisfied 
with the results presented by Kawabata & Mineshige 2009. The most important finding of these similar 
solutions is that, our accreting flow is affected not only by mass-loss but also by loss of energy by the 
wind/outflow. 

There are some limitations in these solutions. One of them is that the accretion flow with conduction 
is a single-temperature structure. So, if one uses a two-temperatures structure for the ions and electrons 
in the discs, it is expected that the ions and electron temperatures decouple in the inner edge, which will 
modify the role of the conduction. The anisotropic character of conduction in the presence of magnetic 
field is the other limitation in our solution (Balbus 2001) 

Although our preliminary similarity solutions are so simplified, they obviously improve our under- 
standing of the physics of hot accretion flow around a compact object. To have a more realistic picture 
of a hot accretion flow, a global solution is needed rather than the self similar solution. In our future 
studies we intend to investigate the effect of thermal conduction on the observational appearance and 
properties of a hot magnetized flow. 
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